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ABSTRACT 

The minimal cooling paradigm for neutron star cooling assumes that enhanced cooling due to neu- 
trino emission from any direct Urea process, due either to nucleons or to exotica such as hyperons, Bose 
condensates, or deconfined quarks, does not occur. This scenario was developed to replace and extend 
the so-called standard cooling scenario to include neutrino emission from the Cooper pair breaking 
and formation processes that occur near the critical temperature for superfluid/superconductor pair- 
ing. Superfluidity is generally expected to exist in the neutron star interior, and Cooper-pair neutrino 
emission processes, which operate through both vector and axial channels, can dominate cooling in 
the minimal model. Neutron stars that have observed temperatures that are too low for their age 
than in the minimal cooling model for any combination of its parameters will imply that enhanced 
cooling is occurring. Previous studies showed that the observed temperatures of young, cooling, iso- 
lated neutron stars with ages between 10^ and 10^ years, with the possible exception of the pulsar 
in the supernova remnant CTA 1, are consistent with predictions of the minimal cooling paradigm as 
long as the neutron ^P2 pairing gap present in the stellar core is of moderate size. 

Recently, it has been found that Cooper-pair neutrino emission from the vector channel is suppressed 
by a large factor, of order 10^"^, compared to the original estimates that violated vector current 
conservation. We show that Cooper-pair neutrino emission remains, nevertheless, an efficient cooling 
mechanism through the axial channel. As a result, the elimination of neutrino emission from Cooper- 
paired nucleons through the vector channel has only minor effects on the long-term cooling of neutron 
stars within the minimal cooling paradigm. We further quantify precisely the effect of the size of 
the neutron ■^P2 gap and demonstrate that consistency between observations and the minimal cooling 
paradigm requires that the critical temperature for this gap covers a range of values between 

T""" ~ 0.2 X 10^ K up to T^^^ ~ 0.5 x 10^ K in the core of the star. This range of values guarantees 
that the Cooper-pair neutrino emission is operating efficiently in stars with ages between 10^ to 
10^ years, leading to the coldest predicted temperatures for young neutron stars. In addition, it is 
required that young neutron stars have heterogenous envelope compositions: some must have light- 
element compositions and others must have heavy-element compositions. Unless these two conditions 
are fulfilled, about half of the observed young cooling neutron stars are inconsistent with the minimal 
cooling paradigm and provide evidence for the existence of enhanced cooling. 
Subject headings: Dense matter — equation of state — neutrinos — stars: neutron 



1. INTRODUCTION 

Neutron stars older than several minutes cool through a combination of emission of neutrinos from their interiors 
and photons from their surfaces. The former dominates the energy losses during the first 10^ years after birth, but 
the latter is responsible for the observed thermal emissions detected from several neutron stars older than hundreds 
of years. For neutron stars of this age and older, the effective surface temperatures are tightly coupled to the interior 
temperatures, and the observability of these stars therefore depends crucially on the overall rate of neutrino emissions 
from within the star. 
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Historically, th eoretical neut r on star cooling models have fallen into two categories, "standard" cooling or enhanced 
cooling (see, e.g. lPethicld[l99l . iPrakashlfTOM Il998l . iPage et al.ll2006l andlPaee II2009D Stai idard cooling implie s that 
no enhanced neutrino emissions from any direct Urea processes (jLattimer et al. 019911 andU rakash et al.lll992l). due 
either to nucleons or to exotica such as hyperons, Bose condensates or deconfined quarks, occurs. iPage et all (j2004f ) 
(hereafter referred to as Paper I) introduced the minimal cooling paradigm to address the question of whether or not 
there exists observational evidence that firmly indicates that enhanced cooling takes place i n some neutron star s. 

Traditionall y, the standard cooling scenario has included only the modified Urea process (|Friman fc Maxwellll979l : 
iTsurutal 119861 ). However, over the last decade i t has been realized that another important cooling mechanism is 
provi ded by Cooper pair breaking and formation (jFlowers. Ruderman fc Sutherlandlll976t iVoskresenskv fc Senatorovl 
Il987f ). termed as the "PBF" process, which occurs in the presence of superconductivity or superfluidity in dense 
matter. It is generally accepted that superfluidity occurs in neutron star matter, although the magnitudes of the 
superfluid gap energies as a function of density are still uncertain at present. Furthermore, for temperatures near 
the associated superfluid critical temperatures, emission from Cooper pairs dominates the neutrino emissivities in 
many cases. Therefore, in the minimal cooling paradigm, neutrino emission is included both from modified Urea 
and bremsstrahlung processes as well as from the Cooper pair breaking and formation processes. Enhanced cooling 
is assumed not to occur. Any neutron star that has an observed surface temperature too low for its age than that 
predicted by the minimal cooling model for any combination of its input ingredients therefore implies that enhanced 
cooling as defined above has occured in that star. 

Using the rates for Cooper-pair pr ocesses established by previous authors (jYakovlev. Kaminker fc Levenfishlll999t 
iKaminker. Haensel fc Yakovlevll 19991 ). Paper I concluded that with the possible exception of the pulsar in 3C58 ^ and 
CTA 1, all cooling neutron stars for which thermal emissions have been detected are marginally consistent with the 
assumption that enhanced cooling is absent, given the combined uncertainties in ages and temperatures or luminosities. 
Importantly, it was concluded that overall consistency with minimal cooling was possible only if the neutron '^P2 gap 
was similar to our model "a" in Paper I, i.e., with critical temperatures of the order of 10^ K. 

This result was not sensitive to the neutron star mass or to the sizes of the n and p ^Sq gaps. In the case that the 
n ''P2 gap was significantly larger than this value, theoretical models yielded too slow cooling to account for the low 
temperatures of about half of the young neutron stars so far observed in thermal emission. It is also apparent from the 
results of Paper I that overall consistency with minimal cooling requires heterogeneity in envelope compositions for 
the young stars: light element compositions for some and heavy element compositions for others. One of the purposes 
of this paper is to show that, besides not being too large, the n ^P2 gap should also be not too small. Unless the n 
■^P2 lies in a narrow range, and there is heterogeneity in young neutron star envelopes, the implication is that about 
half of the observed young neutron stars with thermal emission have some degree of enhanced cooling. 

Of course, even if the conditions on the neutron '^P2 gap and envelope compositions are met, the overall consistency 
of the minimal cooling paradigm with observations does not necessarily mean that direct Urea processes are forbidden. 
Direct Urea processes may in fact still occur in many stars, but th ose neutron stars would quickly grow too cold to be 
detectable from their thermal emissions (Page fc Applegate1 ll992f ). 

Pairing is expected to ap pear at low Fermi momenta {kp ^ 1.5 fm^^) in th e singlet ^Sp state 
(jBohr. Mottelson fc Pinejll958| ) and in the triplet ^P2-'^F2 mixed channel at higher momenta (jTakatsukeJ Il972l ) . In 
both cases, singlet or triplet, Cooper-pair emission can occur through the vector or the axial chann el. The conserva- 
tion of vector cu r rent in dense superfluid neutron matter was first addressed bv iKunku fc Reddvl (|2004D . Recently, 
iLeinson fc Pered ((2Q063) have demonstrated that the vector part of the Cooper-pair emission in a one-component 
system of paired fermions is supp r essed by a factor of ^ (1/20)(i;_f/c)^ ~ 10""^ relative to the original estimates of 
iFlowers. Ru derman fc Sutherland' ('1 976I ). where vp is the velocity of particles at the Fermi surface and c the speed 
of light. iLeinson fc Perez (2006a) have identifed the key reaso n for this large suppression in one-component matter: 
in their original work. IFlowers. Ruderman fc Sutherland! (|197 6l) employed the bare vertex in the vector weak current; 
however, the use of the bare vertex violates vector current conservation. The conservation of vector current is achieved 
by including collective effects which minimally requires (in order to satisfy the Ward identity) that the correction to 
the bare vertex is calculated to the same order of approximation as the quasi-particle propagator is modified by the 
pairing interaction in the system. 

The large suppression in the vector channel has been confirmed by ad ditional works (jSedrakian. Miither fc Schuckl 
120071 iKolomeitsev fc Voskresenskv I l2008t ILeinson I [20081 : ISteiner fc Reddv 2008). If the total Cooper-pair emission 
were suppressed by this factor, significant changes to results of the minimal cooling model might be anticipated. 
However, Cooper-pair emission occurs through both vector and axial channels, and the axial part of the Cooper-pair 
emission, which is proportional to Vp, is only slightly modified. Also, the rate of Cooper pair emission from the triplet 
configuration, being due to axial currents, is largely unaffected by these new results. A cursory examination of the 
results of Paper I suggests additional reasons that the revised Cooper-pair rates might not have a major influence 
on the minimal cooling scenario. One reason is that the specific heat of matter with superfluidity is not affected by 
the existence or absence of Cooper-pair emission. As a result, the cooling associated with ^Sq neutrons, which is 
confined to the crust, has only a transitory effect that is important during the crust's thermal relaxation (i.e., during 
the first few hundred years at most). Another reason is that the cooling associated with the ^Sq protons in the stellar 
core proceeds predominantly through the axial channel. Nevertheless, the large magnitude of the vector suppression 

^ In view of the recent revisions in its age, the luminosity of this object now lies close to the cooling curves of the minimal paradigm. 
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suggests that a quantitative re-analysis of the cooling is in order to confirm these expectations. 

The purposes of this paper are to (i) assess the effects of the suppression of Cooper pair neutrino emission in the 
vector channel on neutron star cooling and to examine their effects on the minimal cooling paradigm, (ii) critically 
evaluate the extent to which current data are compatible with the minimal cooling paradigm, (iii) highlight the crucial 
role of triplet pairing (found in Paper I) in more quantitative terms, and (iv) emphasize that young neutron stars 
cannot have identical envelope compositions and remain compatible with the minimal cooling paradigm. 

In presenting our results below, the key ingredients of the minimal model in Paper I are retained. In doing so, we 
utilize the same (i) equation of state, (ii) superfluid properties of the relevant components, (iii) envelope composition, 
and (iv) stellar mass as in Paper I. We note that in a ne utron star crust, which is a multi-component system, the 
neutron PBF emissivity might be modified by the lattice (jSteiner fc Reddy'2008). However, this effect has not yet 
been computed in detail. For the purposes of this work, we will assume that the modification in the speed of sound is 
small, and, thus, the vjp suppression in the vector channel for homogeneous bulk matter holds in general. The neutrino 
emissivities from the PBF processes will be altered to those calculated recently so that an apposite comparison with 
Paper I can be made. Where appropriate, new inputs with regard to superfluid gaps and revisons in data from 
observations will be utilized and so indicated. 

2. PBF NEUTRINO EMISSIVITIES 

As the temperature nears the critical temperatu re Tc for pairing, new channels for neutrino emission through the con- 
tinuo us formation and breaking of Cooper pairs (iFlowers. Ruderman fc Sutherland! [l976t IVoskresenskv fc Senatorovl 
Il987| ) become operative. 

The PBF emissivity for both neutrons (i = n) and protons (i — p) in the singlet (j — s) and triplet 
(.7 = t) channels (we only consider m j — 0) can be written in the form (jFlowers. Ruderman fc Sutherland! Il976t 
lYakovlev. Kaminker fc Levenfishlll999f ) 

^ ^ ikBTfM.a,,F, [X{T)/T] , (1) 

where Gp is the Fermi weak-interaction constant, m* are effective masses, pF,i are Fermi momenta, T is the temper- 
ature, Afi, is the number of neutrino flavors, j are factors involving vector and axial coupling constants and Fj are 
functions that depend on the temperature dependent pairing gaps Ai(T) and control the efficiency of the PBF process. 
The functional forms of the gaps Ai(T) are discussed in §3.3 of Paper I. Taking Af^, = 3 

Q^3.51 X 10^1 ^ (P^) X rJa,,F, [A,(r)/T] , (2) 

cm'^ s yrrii J \miC J 

where nii are the bare masses , and the control functions Fg and Ft (for mj = 0) are (|Yakovlev. Kaminker fc LevenfishI 
[19991) 

Fs-f i (3) 



/o (1 + e^)' 

Jo [l + e^y 



Ft^ — I I 7^— T2 (4) 



where y = A,(T)/T, z = Jx^ +y^, and / dn represents the angle averaging procedure detailed in 
lYakovlev. Kaminker fc LevenfishI (jl999f ). The form of the control functions Fg and Ft shown in Figure 13 of Pa- 
per I clearly shows that the PBF process turns on when T reaches Tc, increases its efficiency as T decreases, and 

becomes exponentially suppressed when the gap approaches its maximum size A(T = 0) when T ~ 0.2Tc. 

As we wish to examine the extent to which the cooling curves are altered by PBF emissivities that incorporate the 
conservation of the weak vector current, we collect below the relevant results for the factors Uij. 

Vector current not conserv ed: For the singlet and triplet configurations, the factors ai_^ (including both vector and 
axial parts) have the values (|Yakovlev. Kaminker fc LevenfishI 1 19991 : iKaminker. Haensel fc Yakovlevlll999f ) 



an,s = C^T n + CXn PF,n (^1 + ^ j (^) 

^P,s = C^,p + C\ p fp^ {l + (6) 

a„,t = Cy^n + 2C^^„ (7) 

apd — Cy^p + '^C\ .p , (8) 

where pp.i = P fa I "nii and = m*/mi. 

Ve ctor current c onserved: Enforcing the Ward identity suppresses the vector parts of the PBF emissivi- 
ties (iLeinson fc Perez 2006a). The suppress ion in the vector part of the singlet channel is proportional to 



(4/81)u|,j ( Kolomeitsev fc Voskresenskv 1 120081 ). where VF,i = PF^i/m* are the Fermi velocities. The suppression 




Fig. 1. — Rep resentative neutr on ^So pairing critical temperature Tc as a function of the neutron Fermi momentum kp from vari- 
ous calculations: Takatsuka' (11972!) ("T"), a calculation within the BCS approximation; IWambach, Ainsworth & Pines! (119931) ("WAP"), 
[Chen, ct al. (1993) ("CCDK"), and Schwcnk, Friman & B rown (2003) ("SFB"), which take into account medium polarization beyond the 
BCS level; , Gezerlis fc Carlsoa (|200S') ("GC") and Gandolfi. et al.. (.2008) ("GIPSF") denote Quantum Monte-Carlo results (as these latter 
results, marked as triangles and stars, respectively, are restricted to relatively low densities, we have extrapolated them to higher densities). 
In all cases, Tc was obtained from the zero temperature gap A(0) through the standard BCS result as kgTc = 0.57A(0). The vertical 
dotted line shows the location of the crust-core boundary. 

of the vector part in the tri plet channel is also strong, so we set the contribution from the vector part to zero as done 
in iLeinson fc: Perez! (|2006b( ) . Consequently, the quantities modified to 

-P.S=CI., (^) + C%, V\, (l + (10) 

a„,t = 2Ci„ and a^.t = 2C\^. (11) 

The values of the various coupling constants are Cy_„ = 1, Ca,™ = 9A^ Cy.p = 4sin^ 9w — 1, and Ca,p = where 
gA ~ 1-26 and sin^ ~ 0.23. Notice that Cy^ <^ Cy^n, whereas C^ ^ = C'l.n- The fact that pF,p ^ PF.n, imphes 
that neutrino emission from the PBF process involving triplet neutron pairing is intrinsically much more efficient than 

that from singlet proton pairing. 

Although the neutron pairing at high densities involves a mixing of the and ^Fa channels (jTakatsukalll97l . this 
mixing is not taken into account in the triplet PBF emissivities quoted above. In what follows, we will simply refer to 
the pairing of neutrons at high densities as pairing. 

3. EFFECTS OF PAIRING AND PBF EMISSIVITIES 

The most significant revision of PBF neutrino emission is in the case of the neutron ^Sq pairing as emission from 
the '^Pa pairing is suppressed only by about 30%. Comparison of equations (5) through (8) with equations (9) through 
(12) shows that the total emissivity from the proton ^So pairing is essentially unaffected as it is largely dominated by 
the axial channel. Both proton ^Sq and neutron '^Pa pairings occur in the core, which contains more than 90% of the 
star's volume, whereas the neutron ^Sq pairing is essentially restricted to the crust. We now assess how the predicted 
large suppression in the vector channel for the various gaps affects the minimal cooling paradigm. 

3.1. Effect of the neutron ^ Sq gap in the crust 

There are substantial variations in model predictions for the ^Sq gap. For the models we use. Figure [1] summarizes 
predictions for the superfluid critical temperatures Tc as a function of the neutron Fermi momentum kp(n), including 
results of two new calculations performed after Paper I was written. The neutron ^Sq gap is mostly confined to the 
stellar crust, which constitutes only a small fraction of the stellar volume. As a result, effects of the suppression are 
mostly observed during the thermal relaxation of the crust (Lattim er et al. 1 11994: Page 2009). We therefore focus on 
the first 10^ years of evolution during which the effects of suppression are most evident (see Figure [2|). 

We begin with a stellar model in which no pairing is taken into account (curve 1 in Figure [J). In this case, the 
crust's early cooling is driven by neutrino emission from the plasmon process and the neutron-neutron bremsstrahlung 
from the (unpaired) neutrons in the inner crust, with a small contribution from the electron-ion (and a smaller one 
from electron-electron) bremsstrahlung process. 

In the presence of a neutron ^So gap, three effects appear: suppression of the (inner crust) neutron specific heat, 
suppression of the n — n bremsstrahlung, and the onset of the PBF process. We illustrate the influence of these three 
effects in succession in Figure [2j Curve 2 includes only the suppression of the specific heat, which results in a direct 
shortening of the thermal relaxation time of the crust. Curve 3 adds the n — n bremsstrahlung suppression, which 
results in a higher temperature and hence a lengthening of the relaxation phase. Curves 4 and 5 show the total effect 
of the neutron ^Sq pairing by also taking into account the PBF emission process. For comparison, curve 4 is the result 
in which the vector channel contribution is assumed to be unsuppressed whereas curve 5 takes this suppression into 
account. As had been anticipated, the suppression of the vector channel of the PBF process has a significant effect, 
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Log Age [yrs] 

Fig. 2. — Comparison of cooling during the crust thermal relaxation era including various effects from neutron ^So pairing. The ^So gap 
"SFB" from Schwenk. Friman & Brown (2003) has been em ployed. No pairing in the core has been included and the star is a 1.4 Mq star 
built with the EOS of APR ( Akmal. Pandharipande &: Ravenhall.,1998. ) and has a heavy element envelope (see Paper I). Curves 1 through 
5 are explained in the text. 
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Fig. 3. — Early cooling of a neutron star with various neutron ^So gaps. Upper (thick) curves take into account the suppression of the 
vector channel of the PBF process (as in curve 5 of Figure 2) whereas, for comparison, the lower (thin) curves are without the supression 
(as in curve 4 of Figure 2). The dashed curve show the cooling in absence of neutron pairing (as in curve 1 in Figure 2). Notation for gaps 
is the same as in Figure [T] 

but only at early times {t < 1000 years), and results in warmer crusts and increased crust relaxation times. However, 
after crust relaxation is attained, the differences among cooling histories of the various models is very small (this is 
more clearly illustrated below in Figure [8]). 

In order to demonstrate that this result is independent of the specific character of the ^Sq gap, we compare in 
Figure [3] the cooling during the crust thermal relaxation era using the various models for this gap from Figure [TJ 
Although quantitative differences are apparent, the qualitative nature of the effects of suppression are the same for all 
gaps. 

The effects of vector suppression are therefore likely to be important in the interpretation of neut ron star crustal cool- 
ing observ ed in X-r ay transients from accretion-heated neutron stars in Low-Mass X-ray Binaries (ICackett. et al.ll2006l 



2008; Deg enaar. et al. 2009) as the crustal cooling timescale is sensitive to th e neutrino emission rates ()Lattimer et al. I 
1994t iRutledge. et al.ii2002i : IShternin. et al.ll2007t iBrown. fc Cummindl2009f ). Inasmuch as the inclusion of the effects 



of vector suppression increases the cooling timescale, the inferred crust thickness will be overestimated if these effects 
are ignored. Likewise, higher temperatures can be reached in the crust of an accreting neutron star once the vector 
suppression of the PBF rate is taken into account and has important consequence for the triggering of superbursts 
(|Cumming. et al.ll2006( l. 

We therefore conclude that the suppression of the vector channel of the neutron ^So PBF process does not lead to 
a distinguishable effect in the long-term cooling (> 1000 years) of the star. 

3.2. Effects of the neutron ^ P2 and proton ^ Sq gaps in the core 

We now compare the relative efficiencies of the PBF processes from the neutron ■^P2 and/or proton ^So Cooper pairs 
with the modified Urea processes (neutron and proton branches) and nuclcon brcmsstrahlung processes in the core. In 
Figure m we plot the two families of gaps we will consider. As discussed in Paper I, these gaps cover the broad range 
of predicted values from microscopic calculations. Particul arly uncertain is the max i mum size of the n eutron ''P2 gap, 
as well as the range in density over which it is significant (jSchwenk fc Frimanll2004t iBaldo et al.lll99"8l ). 

As an example, we choose the proton ^Sq gap from model "AO" of Figure |4] and we consider three neutron ■^P2 
gaps: one vanishingly small, and the two models "a" and "b" from the same figure, which have maximum values of 
Tc = 10^ K and 3 x 10^ K, respectively. Figure [5] shows the resulting neutrino and photon luminosities in full cooling 
calculations. In the case of a vanishing neutron ■^P2 gap (left panel), the only PBF process occurring in the core 
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Fig. 4. — Critical temperature, Tc, for the proton ^So gaps, left panel, and neutron ^P2 gaps, used in this work. These are the same gaps 
as shown in Figure 9 and 10 of Paper I, respectively, but displayed as functions of the matter density p instead of particle Fermi momenta, 
kp, as in Paper I. Conversion from kp to p was performed with the APR EOS and its corresponding proton fraction. This conversion is 
only weakly dependent on the EOS given the constraints imposed by the conditions of minimal cooling on the EOS. Central de nsiti es of 
stars with masses from 1 to 1.8 Mq are indicated. The grey strip for the neutron ^P2 gaps is the compatibility band discussed in i|3.3l See 
Paper I for references. 




Fig. 5. — Comparison of luminosities from various processes during three realistic cooling histories: photon ("7"), all i/-processes ("Total 
u"), modified Urea and nucleon bremsstrahlung ("MURCA+Brem u"), and PBF ("PBF u") from neutron ■^P2 ("n") and proton ^So 
("p") pairing. PBF neutrino emission from the neutron ^So gap is not shown explicitly as its contribution is always dominated by other 
processes, but is included in the total 1/ luminosity. Suppress ion of the vector channe l of the PBF processes is properly taken into account 
in all cases. In all cases shown, the proton ^So gap is from lAmundsen fc OstgaardI l)1985al ') (model "AO" in Figure |4)l and the neutron 
^So gap from [Schwenk. Friman fc Brownl 1(2003) (model "SFB" in Figure HI. The neutron ''P2 gap is chosen to be vanishingly small 

left panel), from our model "a" (center), or model "b" (right) from Figure |4] The star is a 1.4 Mq star built with the EOS of APR 

Akmal. Pandharipande fc Ravenhalllll998) and has a heavy element envelope (see Paper I). 

is from the proton ^Sq pairing, but due to its intrinsically low efficiency it cannot compete with the modified Urea 
process which is unsuppressed in the inner core where the proton gap vanishes. This model is very similar to the old 
"standard cooling" case. The other two panels in Figure [5] with non- vanishing neutron ^P2 gaps clearly show that the 
PBF processes dominate the total neutrino luminosity as soon as the neutron ■^P2 pairing appears. 

It is worthwhile to note the competition between the proton ^Sq and neutron ^P2 PBF processes which depends on 
the relative sizes of the gaps. For a large neutron ''P2 gap, as our case "b" used to obtain the luminosities in the right 
panel of Figure [SJ the temperature of the entire core drops below Tc in a short time; thereafter, the corresponding 
PBF process is suppressed. When neutrons in the entire core are well into the superfluid phase the PBF process from 

the proton ^Sq gap subsequently drives the cooling, at ages ~ 10'^ yrs, but with a low efficiency. In contrast, when the 
neutron ^P2 model gap "a" is used to obtain the luminosities in the central panel of Figure [S] neutrino emission from 
the neutron PBF process largely dominates the cooling. As noted in Paper I, such gaps as "a" lead to the coldest 
minimal cooling neutron stars. 

3.3. Characterization of the most efficient neutron ^ P2 gaps in the core 

The most efficient pairing configurations, which lead to the coldest neutron stars, a situation explored in Paper I, 
are neutron ■^P2 gaps with Tc values around 10^ K in the largest possible fraction of the core (as in the case of our 
model "a"). In this case, the efficient PBF process from the neutron ^P2 gap dominates the neutrino luminosity at 
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Fig. 6. — A schematic diagram of the neutrino luminosity as a function of temperature for the modified Urea and PBF processes. The 
dotted curve shows the optimal PBF luminosity (i.e., obtainable when a thick enough layer in the core has it temperature close to T^) 
and the dashed curve shows the unsuppressed MUrca luminosity. The values of T™'" and T™^'' are the minimum and maximum values 
of Tc for the neutron gap (which is a function of density) that occur within the star. When the temperature in the core falls below 
jimax ^[^g neutrino luminosity at that point increases to the PBF luminosity, which can be almost two orders of magnitude higher than the 
MUrca luminosity in the optimal case. When the temperature falls further, to below T™'", the neutrino luminosities from both the PBF 
and MUrca process are quenched. (Figure inspired from Figure 20 of Paper I.) 

ages ~ 10" — 10^ yrs (as seen in the central panel of Figure [5]) and results in the coldest young neutron stars within 
the minimal cooling paradigm. 

The schematic illustration in Figure [S] shows the neutrino luminosity as a function of temperature for the modified 
Urea and PBF processes. As long as the temperature is greater than T"^^^, which we define as the maximum value of 
Tc in the core, the modified Urea process drives the cooling. When the temperature falls below T^'^^, the PBF process 
turns on and dominates the cooling, until the temperature drops below T™'", which we define as the minimum value 
of Tc in the core, when both the PBF and modified Urea processes are quenched everywhere in the core. 

The surface temperature at early times is controlled by crustal physics, as described in ij3.11 and is independent of 
the evolution of the core. For the surface temperature to reach the smallest possible values, the value of T"^^^ should 
be large enough for the PBF process to turn on before, or not much later than, the crust isothermalization time. A 
useful reference age is ~ 10'^ years, the estimated age of the youngest observed cooling neutron stars, for the PBF 
process to be fully operating. At later times, if the value of T™'" is too large both the PBF and modified Urea processes 
will turn off before the photon cooling era and the cooling will proceed at a slower pace. This feature is illustrated in 
Figure [7] which shows cooling curves for various neutron ^P2 gaps. The upper solid curve shows cooling for the case in 
which the neutron ■^P2 gap is zero (A = 0). The lower solid curve shows cooling for a neutron ■^P2 gap corresponding 
to case "a", which is close to the most efficient case, giving the lowest temperatures at all ages. The corresponding 
values of T^^^ and T™'" in the 1.4 M© star used in Figure [7] are 10^ K and 2 x 10* K, respectively. The three models 
"0.6a", "0.4a", and "0.2a", in Figure [7] with scaled down gaps show that to obtain the coldest star T^'^^ should be at 
least about 0.5 times the T™'*'' value of model "a", that is, 

T^'^"" > 0.5 X 10^ K. (12) 

With respect to the optimal value of T""", the model "2.0a" in Figure [7] with a scaled up gap shows that the turning 
off of both modified Urea and PBF processes occurs somewhat too early, resulting in warmer stars at ages > 10^ years. 
Consequently, T""" should not be much larger than that of our model "a" for this particular 1.4 Mq star: 

Tc""'" - 0.2 X 10^ K. (13) 

The two bounds in equations [12] and [13] allow us to draw a "compatibility band" for the Tc curve, shown as a horizontal 
grey strip in the right panel of Figure [1] ^ Tc curve will yield the coldest possible minimally cooling neutron star, at 
ages between 10^ to a few times 10'* years, only if it crosses the compatibility band within the density range present in 
the core of the star. Examination of Figure [4] shows that neutron ^P2 gaps resembling those of model "a" provide the 
most efficient cooling for neutron stars of masses close to 1.4 solar masses. 

In order to asses s the role of the gap's density dependence, we also performed cooling simulations using the neutron 
^P2 gap "T" (from fTakat sukal[l 9 72f ) of Figure[4] This gap model is similar to our model "a", but the density dependence 
is different. It has J'™a-x_ q ^ -|^q9 j'mui_Q^ ^^lig minimum being reached at low densities, and thus crosses 

the compatibility band. Cooling simulations utilizing models "a" and "T" result in virtually identical trajectories. In 
contradistinction to the lower bound of equation (jl2p . it is not possible to have an upper bound on T™^^, the only 
condition being that the Tc curve must cross the compatibility band. Notice that a gap with a density dependence 
similar to our models "b" and "c" and a very large T^^^{:^ 0.5 x 10^ K) is also likely to have a very large T""", and 
would violate the compatibility conditions in equations ()12p and (|13p for minimal cooling. This density dependence, 
in par ticular an already large Tc at nuclear matter density, was inspired by the microscopic models of iBaldo et al.l 
()199^ . and the compatibility conditions we find are a strong constraint against such types of gaps. However, in case 
the gap has a vanishingly small Tc a low densities, as the model "T", and grows above 0.5 x 10^ K at higher densities, 
then it will inescapably cross the compatibility band, no matter how large is its T™'^^. It is intriguing that such strong 
constraints on the neutron triplet gap emerge for consistency of the minimal model with data. 
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Fig. 7. — Comparison of predictions of the minimal cooling scenario for variations of the neutron ■^P2 gap. We take as a reference gap 
our model "a" for Figure |4] which we scale by a factor "s"; curves are labelled by this scale factor (the curve 0.0a is the gapless case and 
1.0a represents the standard case "a"). The PBF process ensues when the cooling curve for each gap separates from the upper (gapless) 
trajectory. The standard case 1.0a is seen to repre sent the most effic ient neutrin o emission and the most rapidly cooling case. All models 
are for 1.4 Mq stars built using the EOS of APR l lAkmal. Pandharipande fc Ravenhalliil998i ) with heavy element compositions. 
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Fig. 8. — Summary of observational data on thermal luminosities of isolated cooling neutron stars. Left panel: 12 stars for which a 
thermal spectrum has been clearly detected. Right panel: the PSR in the nebula 3C58 seems to exhibit a thermal component, and could be 
presented in the upper panel, whereas in the case of the other objects a thermal component from the main stellar surface is not detected; 
consequently, the data shown are upper limits on the thermal luminosity. The labels " ?" indicate that a compa ct object has not yet been 
detected in the supernova remnant. Finally, the object RX J0002+6246 in CTBl is possibly not a neutron star (jEsposito et al.ll2008V See 
Appendix A for details. 

4. DISCUSSION AND COMPARISON WITH DATA 

The observational liata against which we will compare our results are displayed in Figure[8]and described in Appendix 
A. An extensive set of full neutron star cooling histories employing the minimal cooling scenario are displayed in 
Figure [5] from which the effects of the vector channel suppression of Cooper-pair emission can be discerned. A total 
of 300 cooling histories are shown and superposed with the data from Figure [H Comparing the left and right column 

panels of Figure [9l it is clear that for times later than the crust thermal relaxation time {t ~ 100 years) the effects of 
suppressing Cooper pair emission from the vector channel are very small. Therefore, the major conclusions that were 
reached in Paper I concerning minimal cooling remain valid. 

In order to highlight the scope and limitations of the minimal cooling paradigm, and to refine its implications, we 
recall its main conclusions as stated in Paper I. Given the uncertainties about the major physical ingredients, namely, 
the sizes of the various pairing gaps and the chemical composition of the envelope, the cooling curves of the minimal 
model appear to be consistent with data of nearly all neutron stars with observed thermal emission, as long as the 
coldest possible trajectories under the minimal cooling paradigm occur in some cases. The compatibility of nearly all 
neutron stars with minimal cooling requires strong conditions on the size of the neutron '^P2 gap and the chemical 
composition of the stellar envelopes, as discussed below. 

A few exceptions exist, notably the pulsar in CTA 1, together with the unobserved emissions from several young 
supernova remnants that might be reasonably expected to contain neutron stars, as they evidently lie below any of the 
minimal cooling histories displayed. These few objects, if confirmed, provide us with the most compelling evidence in 
favor of the occurrence of enhanced neutrino emission beyond the minimal cooling paradigm. 

In the case that the Tc curve crosses the compatibility band as in our models "a" and "T" , agreement with data can 
be achieved if the warmest stars, such as the pulsar in Puppis A and PSR 1E1207-52, have light element envelopes, 
while the coldest ones, like PSR J1119-6127, the pulsars in 3C58 and Vela, and PSR J1124-5916, have heavy element 
envelopes. This agreement is lost, however, if the value of T™™ is too large, i.e., greater than about 0.2 x 10^ K (see 
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Cooling without PBF vector channel suppression Cooling with PBF vector chanel suppression 




Log t [yrs] Log t [yrs] 



Fig. 9. — Comparison of predictions of the minimal cooling scenario with data; all models are for 1.4 Mq stars built using the EOS of 
APR (Akmal, Paridharipandc & Ravcnhall 1998). In the right panels the suppression of the vector channel in the Cooper-pair neutrino 
emission is fully taken into account whereas, for comparison, in the left panels the supression has been omitted. In each row, the two panels 
have the same neutron gap, from a vanishing gap in the upper row to our model gaps "a" and "b" (following the notations of Figure 
10 in Paper I) in the next two rows. In each panel two sets of cooling trajectories, either with light or with heavy element envelopes, are 
shown which include 25 curves corresponding to 5 choices of the neutron ^Sq and of the proton ^Sq gaps covering the range of predictions 
about the sizes of these gaps. 

equation [13]), as in our models "b" and "c" . In the extreme case that the neutron gap is vanishingly small and also 
that all observed young cooling neutron stars have light element envelopes, then nearly all of them, with the possible 
exception of PSR B0538+2817, are observed to be too cold to be compatible with minimal coohng predictions. In 
the less extreme possibility of a heterogeneity in chemical composition and a vanishingly small neutron gap, we 
still find that more than half (seven out of twelve) of the observed young cooling neutron stars are too cold to be 
compatible with minimal cooling. (Notice that among the remaining five, out of twelve stars, the compact objects in 
Gas A and the Grab still have only upper limits.) If these conditions on the Tc curve are not satisfied for a particular 
model of superfluidity in dense matter, then that model also requires enhanced cooling beyond the minimal cooling 
paradigm. These results highlight the importance of the n ■^P2 gap in more precise terms than discussed in Paper I. 

Our conclu sion regarding the need f or heterogenity in the chemical composition of the atmosphere is consistent with 
the results of lKaminker. et al. I ()2006D . who had to employ both light and heavy element atmospheres in their cooling 
models to match the data of most stars. 

That it is apparently possible to explain the majority of thermally-emitting neutron stars with the minimal cooling 
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paradigm led iKlahn. et al.l (|2006D to propose the so-called direct Urea constraint which would effectively place limits 
on the density dependence of the nuclear symmetry energy. This follows from the fact that higher mass stars have 
larger central densities, so that the few stars that appear to be too cold could have larger than average neutron 
star masses. The nearly two d ozen well-measured neutron star masses are mostly concentrated in the range 1.25-1.5 
Mq (Latti mer fc Pr akash 2005), although there are indications that some neutron stars might extend be yond the upper 
end o f this range. Mass-depende nt cooling in the context o f dir ect Urea processes was sug gested by iLattimer et al.l 
(|1991h . and has been exploited bv iPage fc Applegate I ()1992[ ) and lYakovlev fc Pethickl ()2004f ) to account for the gross 



features of cooling data. If rapid cooling is, in fact, due to the onset of the direct Urea process involving nucleons, 
as opposed to hyperons, kaon condensates, or deconfined quarks, a constraint on the density dependence of the 
nuclear symmetry energy could therefore result. The nuclear symmetry energy would have to grow relatively slowly 
with density so that the threshold density for the direct Urea process to occur exceeds the central density of most 
neutron stars, few of which could be very massive. The inferred constr aints on the sy mmetry energy may also depend 
significantly on the presence of quartic terms in the symmetry energy (jSteiner II2006I ). 

However, if any of the assumptions in the above logical sequence are violated, the direct Urea constraint becomes 
invalid. Envelope compositions and the neutron ■^P2 gap might be outside the required ranges. Perhaps even in massive 
neutron stars cooling due to the direct Urea process is abated due to superfluid-induced suppression. Furthermore, 
in the exceptional cases in which rapid cooling is suspected, it will be difficult to prove that the nucleon direct Urea 
is responsible as all proposed rapid cooling emission rates are orders of magnitude greater than modified Urea or 
Cooper-pair emission rates. It should also be noted that there are at present no reliable mass measurements of any 
thermally-detected cooling neutron stars, so that mass information deduced from binary pulsars might not be relevant. 

5. SUMMARY AND CONCLUSIONS 

Neutrinos emitted in the continual Cooper-pair breaking and formation (PBF) processes that occur near the critical 
temperatures of superfluid neutrons and superconducting protons in neutron star interiors are an integral part of the 
minimal cooling paradigm (Paper I). These processes were not generally included in the traditional standard cooling 
scenario, but are nevertheless expected to occur whether or not rapid cooling occurs. If all neutron star cooling data 
are found to be consistent with the predictions of the minimal cooling model, (in which rapid cooling by direct Urea 
processes occurring through reactions involving nucleons, hyperons, quarks, or Bose Condensates were eliminated by 
design) then there is no reason to invoke rapid cooling. Either rapid cooling does not occur or it is suppressed by 
superffuidity. In either case, co nstraints about th e propert ies of dense matter en sue. 

However, beginning with iLeinson fc Pered 1200631). several auth ors (jSedrakian. Miither fc Schuckl 120071 : 
iKolomeitsev fc Voskresenskv I l2008t iLeinson I l2008t ISteiner fc Reddvl |2008[) have recently established that the PBF 
neutrino emissivity from the vec t or cha nnel is suppressed by a factor of ~ 10"'^ from the original estimates of 
iFlowers. Ruderman fc Sutherland ()1976l ). who overlooked the conservation of vector current in weak interactions. 
In view of this development, we have performed extensive calculations to assess the effects of the severely reduced 
rates in the vector channel on the long-term cooling of neutron stars incorporating the revised PBF rates. 

Our analysis leads us to conclude that the long-term cooling of neutron stars as envisaged in the minimal cooling 
paradigm is not significantly affected by the proposed large reductions in the vector channel of the PBF neutrino 
emissivities. The reason is that Cooper-pair emission occurs through both vector and axial channels, and the axial 
part of the Cooper-pair emission is barely affected. The axial channel of the PBF emissivites controls the long-term 
cooling, whereas the vector channel emissivities are important for the cooling behavior of neutron stars with ages 
less than a few hundreds of years for which observational data are not yet available. With the exception of few 
candidates (e.g., CTA 1 and perhaps unobserved stars in other remnants), the minimal cooling paradigm is consistent 
with observations of neutron star temperatures and ages, but, as noted in Paper I, only if certain combinations of 
light- or heavy-element envelopes and sizes of neutron triplet pairing gaps in the star's core exist. 

For the minimal cooling model to be consistent with data, it is necessary for the most efficient cooling possible to 
occur. Therefore, our analysis places a stringent requirement on the critial temperature for neutron superfluid pairing 
in the triplet channel in the core of the star. It is required that the Tc curve crosses the compatibility band described 
in Figure m and therefore the triplet gap must both become larger than Tc — 0.5 x 10^ K at some moderately large 
density (but less than the central density of the neutron star), and also become smaller than Tc — 0.2 x 10^ K (either 
at high density, i.e., in the center of the star, or at low density, i.e. in the outermost part of the core). If the gap does 
not fulfill these conditions, then observations point to the occurence of enhanced cooling in at least half of the young 
isolated cooling neutron stars observed so far. 
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APPENDIX 
OBSERVATIONAL DATA 

The observational data plotted in Figures [5] and [5] in the text deserve some comments. We have considered the same 
objects and data as reported in Paper I, but with modifications for the first five objects discussed below. We have also 
included some additional objects (see the last four objects below) for which data has become available since Paper I 
was published. Salient aspects of the modifications and new data are summarized here. 

In many cases, the only way in which the ages of the objects can be estimated is by their spin-down ages 2P/ P. 
However, there are many stars for which kinematic ages have also been obtained. Comparing kinematic and spin-down 
ages for the same objects often reveals discrepancies of factors of three or more in both directions. For example, SNR 
N157B has a determined age of less than 2000 years but a pulsar spin-down time of 5000 years, and PSR B1757-24 
has an estimated age less than 39000 years and a spin-down age of 16000 years. For this reason, we have attached 
nominal error bars of a factor of three to those objects for which we only have spin-down information. To ignore this 
uncertainty could be misleading. 

PSR RX J0205+6449 in the nebula 3C58 

The supernova remnant 3C58 has been associated with the historical supernova SN 1181 (iSteohenson fc GreenlfToMl 
providing us with a pulsar age of 828 years. However, this proposed association has been challenged in recent years. 
Measurements of the remnant's expansion velocity compared to its size (at an assumed distance of 3.2 kpc) require an 
older origin. Optical observations of filament expansion (angular expansion) and filament line width (radial expansion) 
imply an age of 3000-4000 yrs and 2700-3500 yrs, respectively (Fcscn ct al. 2008), whereas radio measurements of the 
filament expansion would imply an age of 4500-7000 yrs (Bictenholz 2006). More over, model ing of the energetics of 
the pulsar wind nebula and its expan sion need an age of 2000-3000 yrs (Chevalieil[2004 l2005f ) or even 2000-5000 yrs 
from the recent X-ray observations of lGotthelf. Helfand. fc Newburghl (|2007h . These ages are more in accord with the 
pulsar spin-down age of 5400 yrs. We will, henceforth, adopt a conservative age of 2000-7000 yrs, keeping in mind the 
possible SN 1181 association and the corresponding age of 828 yrs. 

PSR B 0833-45 in the SNR Vela 

The kinematic age of the SNR we adopted in Paper I was taken from lAschenbach. Egger fc Triimpeil ()1995l ) who 
employed a distanc e of 500 pc (the PS R dispersion measure distance). Correcting this dista nce to 300 pes ( the V LBI 
parallax distance of iDodson et al.ll200i!H gives us an age of 5400-16000 yrs as emphasized bv iTsuruta et al.l ()2009i) . 

RX J0007. 0+ 7302 in the SNR CTAl 

We notice that t he Fermi Gamma- Ray space telescope has detected pulsations at 315 ms from this compact source 
(jAbdo et al.l l2008l ) and measured its spin-down age as 1.4 x 10'' yrs, in agreement with the estimated age of the 
remnant, 1.4 — 2.2 x 10* yrs. 

PSR 1055-52 

This pulsar's spin-down age is 10^'^^ yrs instead of 10^'*'^ yrs which was erroneously adopted in Paper I. 

PSR 1E1207-52 

A lower limit on the spin-down age of this pulsar has been found to be 27 Myrs (jGotthelf fc Halpernll2007l ) while the 
associated supernova remnant kinematic age is about 7000 yrs (within a factor 3). As in Paper I, we use the kinematic 
age from the SNR. 

PSR B2334-t-61 in the SNR G 114.3 -1-0.3 
Surface thermal emission has been detected in an XMM-Ncwton observation (McGowan et al."2006'). The distance 
estimate from the pulsar dispersion measure is 3.1^^ q kpc, but HI absorption (Yar-Uy aniker e t al. 2994) indicates a 
much shorter distance of less than 1 kpc. Fits of the spectrum with a magnetized hydrogen neutron star atmosphere 
(NSA) gives T^^^ = 0.65 x 10^ K and R^^^ ~ 13 km, whereas a blackbody (BB) fit gives T^^ = 1.6 x 10^ K and 
i?^^ ~ 1.6 km. In both distance of 3.1 kpc was assumed. The hydrogen atmosphere model seems to give better 

fits, but the BB fit would be preferred in case of the smaller distance. However, the uncertainty on the luminosity is 
dominated by the uncertainty in the distance and with the above values, we obtain Loo = 0.5 x 10^^ — 5 x 10'^^ erg 
s^^ using the NSA values, whereas the BB values give an Loo range included within the NSA results. The age of the 
SNR, a few times lO'' yrs, is uncertain as is the association with the PSR whose spin-down age is 4 x 10* yrs. We will 
adopt an age of (^ — 3) x tsd, i-e., with the same uncertainty factor as we did in Paper I. 



12 



PAGE, LATTIMER, PRAKASH, & STEINER 



PSR 31119-6127 

This high- field pulsar shows a thermal speetrum which when fitted by a blackbody spectrum gives = (1.35 — 
2.25) X 10^3 erg 8"^ NSA fits yield L'^^^ = (1.05 - 2.05) x 10^3 erg s'^ (|Safi-Harb &: Kumaj[200l . The NSA fit is 
compatible with emission from the entire neutron star surface whereas the BB fit implies a much smaller radius (~ 3.5 
km) at an estimated distance of 8.5 kpc. However, the distance is uncertain and the exceptionally high pulsed fraction 
(~ 75 %) (IGonzalez et al.l[2005l) may favor a small thermally emitting region. Given this and because and 
are similar, we adopt = (1.05 — 2.25) x l O^^ erg s~^ for the t hermal luminosity. The spin-down age is tgd = 1600 
yrs. As the measured braking index n — 2.9 (jCamilo et al.|[2Q00l ) is compatible with magneto-dipolar breaking, and 
the characteristics of the a ssociated SNR and the Pulsar Wind Nebula (PWN) are compatible with an age close to tgd 
(jSafi-Harb &: Kuniai1l2008l ). we adopt an age estimate of (0.5 — 2) x tsd- 

PSR BOSS 1+21 in the Crab nebula 

Phase resolved X-ray spectroscopy (jWeisskopf et al.ll2004l ) gives the most stringent upper limit on the surface thermal 
luminosity, L^o < 2.8 x 10^* erg s~^. However, there is no evidence of a thermal component in the spectrum which is 
completely dominated by magnetospheric emission. 

SNRs G065.3+5.7 and G074.0-8.5 

In their extensive search for compact objects in nearby SNRs, I Kaplan et al.l ()2006D provide six new upper limits on 
the thermal luminosity of possibly existing (but undetected) compact stars. Four of these limits are quite high so we 
will retain only the two lowest ones: Loo < 10^^ erg s~^ for SNR G065.3-f 5.7 (Cygnus Loop) and L^o < 0.7 x 10^^ erg 
s^^ for G074.0-8.5. The estimated ages of the remnants are (2 — 2.5) x lO"' yrs and (0.6 — 1.2) x 10** yrs, respectively. 
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